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Abstract 

The (complex) two-qubit systems comprise a 15-dimensional convex set and the real two-qubit 
systems, a 9-dimensional convex set. While formulas for the Hilbert-Schmidt volumes of these 
two sets are known — owing to recent important work of Sommers and Zyczkowski (J. Phys. 
A 36, 10115 [2003]) — formulas have not been so far obtained for the volumes of the separable 
subsets. We reduce these two problems to the determination of certain functions (freai(u) and 



fcompiex(u)) °f a single variable u = y P22P33 ' wnere P 1S ^ e corresponding 4x4 density matrix, 
and the pa's (i = 1,...4) its diagonal entries. The desired separable volumes are, then, of the 
form V sep / real = 2 jac rea i(u)f rea i(u)du and V comp i ex/real = 2 jac comp i ex (u)f comp i ex {u)du. Here 
jac denotes the corresponding (known) jacobian function. We provide estimates of the two sets of 
/0*)'s and Vs. 
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I. INTRODUCTION 

In a pair of major, skillful papers, making use of the theory of random matrices Q], 
Sommers and Zyczkowski were able to derive explicit formulas for the volumes occupied by 
the d = (n 2 — l)-dimensional convex set of n x n (complex) density matrices (as well as the 
d = — 1 ^ TO+2 ) -dimensional convex set of real (symmetric) n x n density matrices), both in 
terms of the Hilbert-Schmidt (HS) metric |2j - inducing the flat, Euclidean geometry - 
and the Bures metric |j| (cf. Q]). Of course, it would be of obvious considerable quantum- 
information-theoretic interest in the cases that n is a composite number, to also obtain 
HS and Bures volume formulas restricted to those states that are separable — the sum 
of product states — in terms of some factorization of n j^J. Then, by taking ratios - 
employing these Sommers-Zyczkowski results — one would obtain corresponding separability 
probabilities. (In an asymptotic regime, in which the dimension of the state space grows to 
infinity, Aubrun and Szarek recently concluded [f| that for qubits and larger-dimensional 
particles, the proportion of the states that are separable is superexponentially small in the 
dimension of the set.) 

In particular, again for the 15-dimensional complex case, n = 4 = 2 x 2, numerical 
evidence has been adduced that the Bures volume of separable states is (quite elegantly) 
2-15(3^1) ~ 4.2136 • 10" 6 0, Table VI] and the HS volume (5V3)~ 7 w 2.73707- 10" 7 eq. 
(41)]. Then, taking ratios (using the corresponding Sommers-Zyczkowski results), we have 
the derived conjectures that the Bures separability probability is 1680 ^~ 1 ) ~ 0.0733389 and 
the HS one, considerably larger, ^'^V 13 ^ ~ 0.242379 ^, eq. (43), but misprinted as 5 3 
not 5 4 there]. (Szarek, Bengtsson and Zyczkowski — motivated by the numerical findings 
of 0, 0| have recently formally demonstrated "that the probability to find a random 
state to be separable equals 2 times the probability to find a random boundary state to be 
separable, provided the random states are generated uniformly with respect to the Hilbert- 
Schmidt (Euclidean) distance. An analogous property holds for the set of positive-partial- 
transpose states for an arbitrary bipartite system" [lCJ] (cf. Hj). These authors also noted 
jiol p. L125] that "one could try to obtain similar results for a general class of multipartite 
systems" . In this latter vein, recent numerical analyses of ours give some [but certainly not 
yet conclusive] indication that for the three-qubit inseparable states, there is an analogous 
probability ratio of 6 — rather than 2.) 
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However, the analytical derivation of (coniecturally) exact formulas for these HS and 

n nn 

Bures (as well as other, such as the Kubo-Mori |12j and Wigner-Yanase |8j, Il3j ) separable 
volumes has seemed quite remote — the only analytic progress to report so far being cer- 
tain exact formulas when the number of dimensions of the 15-dimensional space of 4 x 4 
density matrices has been severely curtailed (nullifying or holding constant most of the 15 
parameters) to d < 3 QQ (cf. H). Most notably, in this research direction, in Q, 
Fig. 11], we were able to find a highly interesting/intricate (one-dimensional) continuum 
(— oo < j3 < oo) of two-dimensional (the associated parameters being 6 1( the mean, and 
<7g, the variance of the Bell-CHSH observable) HS separability probabilities, in which the 
golden ratio jl7[ was featured, among other items. (The associated HS volume element - 
1 dfidbqda^ — is independent of b± and a 1 in this three-dimensional scenario.) Further, 
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14j, building upon work of Jakobczyk and Siennicki 18[, we obtained a remarkably wide- 



ranging variety of exact HS separability (n = 4, 6) and PPT (positive partial transpose) 
(n = 8, 9, 10) probabilities based on two-dimensional sections of sets of (generalized) Bloch 
vectors corresponding tonxn density matrices. 

The full d = 9 and/or d = 15, n = 4 real and complex two-qubit scenarios are quite 
daunting — due to the numerous separability constraints at work, some being active [binding] 
in certain regions and in complementary regions, inactive [nonbinding] . "The geometry of 
the 15-dimensional set of separable states of two qubits is not easy to describe" jlOj, p. L125]. 
We will seek to make substantial progress in these directions here, by recasting both these 
problems within a one-dimensional framework. 

To proceed in our study, we employ the (quite simple) form of parameterization of the 



density matrices put forth by Bloore 



19 . |20j some thirt y y ears 



a number of other possible parametrizations 21 
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Of course, there are 
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271 ] . a number of which 



we have also utilized in various studies |28j, |29| to estimate volumes of separable states. Our 
greatest progress at this stage, in terms of increasing dimensionality, has been achieved with 
the Bloore parameterization — due to a certain computationally attractive feature of it, 
allowing us to decouple diagonal and non-diagonal parameters — as detailed shortly below.) 
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II. BLOORE PARAMETERIZATION OF THE DENSITY MATRICES 

The main presentation of Bloore was made in terms of the 3 x 3 (n = 3) density 
matrices. It is clearly easily extendible to cases n > 3. The fundamental idea is to scale the 
off-diagonal elements (pij,i ^ j) of the density matrix in terms of the square roots of the 
diagonal entries (pa). That is, we set (introducing the new [Bloore] variables Zij), 

Pij — \/ PiiPjjZij- (1) 

This allows the determinant of p (and analogously all its principal minors) to be expressible 
as the product (\p\ = AB) of two factors, one (A = U.f =1 pu) of which is itself simply the 
product of (nonnegative) diagonal entries (pa). In the real n = 4 case under investigation 
here — we have 

B — (^34 — l) z l2 + ^ (Z\4 (Z24 ~ ^23^34) + z 13 ( z 23 ~~ ^24^34)) z 12 ~ z 23 ~~ ^24 — Z 34+ (2) 
Z 14 (^23 ~~ l) + Z 13 { Z 24 ~ l) + 2223^24^34 + 2 Z 13 Z 14, ( z 34 ~ z 23 z 2l) + 1, 

involving (only) the z^s (i > j), where Zji = |3, eqs. (15), (17)]. Since, clearly, the 
factor A is positive in all nondegenerate cases (pa > 0), one can — by only analyzing B 
- essentially ignore the diagonal entries, and thus reduce by (n — 1) the dimensionality of 
the problem of finding nonnegativity conditions to impose on p. This is the feature we will 
seek to maximally exploit here. A fully analogous decoupling property holds in the complex 
case. 

It is, of course, necessary and sufficient for p to serve as a density matrix (that is, an 
Hermitian, nonnegative definite, trace one matrix) that all its principal minors be nonneg- 
ative [3^. The condition — quite natural in the Bloore parameterization — that all the 
principal 2x2 minors be nonnegative requires simply that —1 < z%j < l,i ^ j. The joint 
conditions that all the principal minors be nonnegative are not as readily apparent. But for 
the 9-dimensional real case n = 4 — that is, Q{pij) = — we have been able to obtain one 
such set, using the Mathematica implementation of the cylindrical algorithm decomposition 
j^j]. (The set of solutions of any system of real algebraic eq uations and inequalities can 
be decomposed into a finite number of "cylindrical" parts |32j.) Applying it, we were able 
to express the conditions that an arbitrary 9-dimensional 4x4 real density matrix p must 
fulfill. These took the form, 212, 213, Z14 G [—1,1] and 

^23 ^ [Z 23 , ^23], ^24 £ [-^24; -^24]; Z 34 £ [-^34)^34]) (3) 



where 



Z 23 = z 12 z 13 ± y 1 - 4,y 1 - zf 3 , Z 24 = 212214 ± y 1 - 2? 2 V 1 - «14> ( 4 ) 

7 ± _ z 13 z 14 — z 12 z uZ23 — Z 12 Z 13 Z 2 4 : + 223^24 ± S 

34 ~~ 1 _ ^2 ) 

1 z 12 



and 



S — y-1 + 2i2 + 2^3 — 2Zi22l3^23 + Z 23y ~^ + ^12 + ^14 ~~ 22i22l4^24 + ^24- (5) 

Making use of these results, we were able to confirm via exact symbolic integrations, the 
(formally demonstrated) result of Zyczkowski and Sommers [2] that the HS volume of the 
real two-qubit (n = 4) states is ~ 0.0016106. (This result was also achievable through 
a somewhat different Mathematica computation, using the implicit integration feature first 
introduced in version 5.1. That is, the only integration limits employed were that z^ 6 
[— 1, 7^ j — broader than those in (jSJ) — while the Boolean constraints were imposed 
that the determinant of p and one [all that is needed to ensure nonnegativity] of its principal 
3x3 minors be nonnegative.) 



A. Determinant of the Partial Transpose 



owever, when we tried to combine these integration limits (JHJ) with the ( Peres- Horodecki 



34, 



35| n = 4) separability constraint that the determinant (C = |ppt|) of the partial 
transpose of p be nonnegative |36l Thm. 5], we exceeded the memory availabilities of our 
workstations. In general, the term C — unlike the earlier term B — unavoidably involves 
the diagonal entries (pa), so the dimension of the accompanying integration problems must 
increase, it would seem, we initially thought — in the 9-dimensional real n — 4 case from 6 
to 9. 

However, we then noted that, in fact, the dimensionality of the required integrations must 
only essentially be increased by one (rather than three), since C turns out to be (aside from 
the necessarily nonnegative factor of A) expressible solely in terms of the (six, in the real 
case) distinct z^'s and the square root ratio 



^ _ pupu ^ 
P22P33 
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(Considering [i as fixed, this is the equation of an hyperboloid of one sheet [37J, p. 227].) 
That is, 

C= \p PT \ = A^-zl^ 4 + 2z M (z 12 z r3 + V 3 + s^ 2 + 2z 2 3(z 12 Z24 + z 13 z u ) p-z^, (7) 
where 

S = { Z 34, — l) Z 12~ 2 (^14^23 + Z13Z24) Z34Z12 — zl 3 + zf 4 zl 3 + {z\% — l) 2^4 — ^34 ~ 2Zi3Zi42;23-2 ; 24 + 1 • 

C is, thus, a quartic/biquadratic polynomial in terms of ji (cf. (23, Q). (Clearly, the 
difficulty of the two-qubit separable volume problem under study here is strongly tied to 



the high [fourth] degree o: 
can be reduced to 2 (cf. 



C in /i. By setting either z 14 = or z 23 = 0, the degree of C 
) In the complex case, C once again assumes the form of a 
quartic polynomial in /jl. So one encounters, in that setting, thirteen-dimensional integration 
problems rather than fifteen-dimensional ones. 



III. ANALYSES 



So, the problem of determining the separable volumes can be seen to hinge on (in the real 
case), a seven-fold integration involving the six (independent) ^'s and \i. However, such 
requisite integrations, allowing /i to vary (or even holding \i constant at various values, thus, 
reducing to six-fold integrations), did not appear to be exactly/symbolically performable 
(using version 5.2 of Mathematica). 



A. Estimation of the univariate functions f r eai(p) and f complexity) 

Thus, to make further progress, it seemed necessary, at this stage, to employ numeri- 
cal methods (not excluding the possibility that exact solutions might, at some point, be 
revealed) . 

We proceeded along two parallel courses, one for the 9-dimensional real two-qubit case 
and the other for the 15-dimensional complex case. We sought those functions freai(^) and 
f compiexif 1 ) that would result from imposing the conditions that the expressions A, B and C 
(as well as a principal 3x3 minor of p), along with their complex counterpart expressions, 
be simultaneously nonnegative. (The satisfaction of these joint conditions ensures that we 
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are dealing precisely with separable 4x4 density matrices.) It was evident that the relation 
f(fx) = /(-) must hold, so we only studied the range /i G [0, 1]. Dividing this unit interval 
into 2,000 equal nonoverlapping subintervals of length ea ch, we sought to estimate the 
/(/i)'s at the 2,001 end points of these subintervals. 

This required (/i being fixed at these end points) numerical integrations in 6 and 12 
dimensions. For this purpose, we utilized the Tezuka-Faure (TF) quasi-Monte Carlo pro- 



4£j, we had extensively used in our earlier studies of separability probabilities 



cedure 

[3, sl- For each of the 2,001 discrete, equally-spaced values of \i we employed the same set 
of 37,000,000 Tezuka-Faure six-dimensional points in the real case and, similarly, the same 
set of 25,000,000 twelve- dimensional points in the complex case. (The Tezuka-Faure points 
are defined over unit hypercubes [0, l] n , so in our computations, we transform the Bloore 
variables accordingly. We plan to continue to add such points to our [real and complex] 
analyses.) 

In Figs. Hand El we show the results of this procedure. There were some slight deviations 
from monotonicity 44] (presumably due to limited sample sizes) in the vicinity of \i = 1 for 
both functions. 

In the real case, our estimate of known Hilbert-Schmidt volume of (separable plus non- 
separable) states [2j, ~ 0.0016106 was larger by only a factor 1.00006. So, we would 
expect our companion estimates of f r eai(f>), at each of the 2,001 sampled points, to be 
roughly equally precise. (Let us note that / rea z(0) = /complexity = 0-) I* 1 the complex case, 
our estimate of the known 15- dimensional volume, OK1Q 7r K 'Lnn ~ 1.12925 ■ 10~ 6 was smaller 

' 851350500 

only by a factor of 0.99965. (As instances of specific values, based on independent analy- 
ses using still larger numbers of TF-points, we obtained estimates of / reo j(l) = ^fff§||p ~ 
114.62368, f real {\) = « 74.108728, f r U°Au) = Hif « 88.312344, f complex (l) = 

387.33307366, f complex {\) = 180.6046580, f complex (a Au ) = 251.157815860, where a Au = 
denotes the golden ratio |l2j. Exact characterizations of f re ai(fj) and f C om P iex(^) would, of 
course, be of great interest, in particular, for the possibility that they might yield exact 
volume results.) 

To estimate the desired separable volumes (V sep ) themselves, one must perform the cal- 
culations, 

Vsep/real = 2 / j ac realfreal{ij)d^ (8) 
JO 
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FIG. 1: Estimation of f r eai(n) based on the third-order interpolation of 2,001 points (/i), the value 
at each such point being based on six-fold numerical integrations employing the same (for each fi) 
set of thirty-seven million Tezuka-Faure points 

f complex t ) 
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FIG. 2: Estimation of f complexity) based on the third-order interpolation of 2,001 points the 
value at each such point being based on twelve-iold numerical integrations employing the same (for 
each ji) set of twenty-five million Tezuka-Faure points 



and 



V, 



sep j complex 



j dC-complex f complex ( H> ) dfl . 



(9) 



B. Jacobians for the transformations 



Now (Fig. EJ), 

// (12 ((/i 2 + 2) Gu 4 + IV + 8) /i 2 + 1) log(n) - 5 (5/i 8 + 32/i 6 - 32^x 2 - 5)) 



jaC rea l{n) = 

and (Fig. EJ) 



1890 (n 2 - 1) £ 



(10) 
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FIG. 3: Plot of the jacobian function jac rea i(fj,), given by lfTU)l 
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FIG. 4: Plot of the jacobian function jac comp / e3; (^), given by (fTTj) 



where 



j dCcomplex (/^) 



1801800 O 2 



-V, 



(11) 



y = 363y^ ±4 + + 48363/i iU + 42875// - 42875/Z 5 - 48363/x 4 - 9947/T - 363 

-140 (fi 14 + 49// 2 + 441/x 10 + 1225/x 8 + 1225/i 6 + 441/z 4 + 49/i 2 + l) log(//). 
We have that 



jac rea i({i)d{i 



71 



and 



2293760 
1 



4.30281 • 10~ 6 , 



ja-Ccomviexi^du = ~ 4.95536 • 10 

o J com P lex ^ J p 2018016000 



-10 



(12) 



(13) 



(The smallest value of \i for which jac rea i(fj,) = that we were able to find was 0.9685588023, 
while for iac comp i ex {jj) we found 0.8395384257.) 
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We obtained the jacobian functions jac reo i{p) and jac comp i ex (p), given in (fTUj) and (JUJ), 
by transformations of, say, P33 to the p variable (and subsequent two-fold exact integrations 
over pu and P22) of the original (three-dimensional) jacobians, involving the diagonal entries, 
for the Bloore parameterizations. These original jacobians were of the form (IT* =1 pjj) fc with 
k = I in the real case, and k = 3, in the complex case. (Of course, by the unit trace 
condition, we must have p^ — 1 — pn — P22 — P33-) 

The direct high-accuracy computation of the desired separable volume integrals (jHJ and 
(jHJ) proves challenging due to the highly oscillatory nature of jac rea i(p) and jac comp i ex (p) 
(given by (jlOj) and (jll)) ) in the vicinity of /i = 1, as indicated in Figs. El and EJ (It might be 
appropriate to sample more points in the vicinity of p. = 1 than in other less problematical 
regions. We have also attempted — without significant success so far — to evaluate these 



integrals using repeated integration by parts 
repeated exact integrations.) 



42j, since the two jacobians in question admit 



C. Volume integrals over p G [0, .95] 

Replacing the upper integration limit of 1 in the integral (jHJ) by .95, we obtained — using 
high precision arithmetic — a result of 0.0006707668 and consequent lower bound on the 
probability of separability of the real two-qubit systems of 0.41647013. (The direct use of 
upper integration limits greater than .95 appeared to lead to unstable results.) For similar 
reasons, replacing the upper integration limit of 1 in the integral © also by .95, we obtained 
a result of 2.327058044 ■ 10~ 7 and consequent lower bound on the probability of separability 
of the complex two-qubit systems of 0.2060707612. 



D. Volume integrals over p G [.95, 1] 

In the immediately preceding analysis, we used upper limits of .95 rather than 1 in the 
integrals (jHJ) and (jH)). To estimate the integrals in the remaining range [.95, 1], we replaced 
the jacobian functions jac rea i(p) and jac comp i ex (p), given in (flUl) and (JUJ), by their 100-degree 
power series expansions about p = 1. (When plotted over [.95,1], both these replacement 
functions gave the appearances of simple downward- sloping lines (Figs. El and EJ).) 

Proceeding in such a manner, again using high-precision arithmetic and summing the 
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approx . 




FIG. 6: 100-degree power series approximation to jac COTO?) / ea ;(/x) about u = 1 

results over the two sets of intervals, we arrived at our final (subject to the availability 
of additional Tezuka-Faure points) estimates V sep / rea i ~ 0.0007298112 and V sep / comp i ex « 
2.625622678 • 10~ 7 . Then, we have prob sep / rea i « 0.45313001 and prob comp i ex / rea i ~ 
0.23250991. (When we compared these several results, based on interpolation — to es- 
timate freaii^) and fcompiexin) (Figs. [I] and |2J) — using third-degree polynomials with those 
using sixth-degree polynomials, we obtained essentially the same set of results.) 

As noted in the introductory section, we had previously hypothesized that V sep / com ni ex = 
(5V3)- 7 « 2.73707- 10" 7 eq. (41)] and prob sep/complex = ^^ii-ia^ w .242379 Q, eq. 
(43), but misprinted as 5 3 not 5 4 there]. The analysis there was based on a considerably 
larger number — 400,000,000 — of Tezuka-Faure points than here. But each point there 
was employed only once for the Peres-Horodecki separability test, while each point here is 
used in 2,000 such tests (with \x ranging over [0,1]). (Relatedly, we had initially suspected 
that if we started checking the Peres-Horodecki criterion for successively larger values of 
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holding the set of z^s given by a Tezuka-Faure point fixed, then if we reached one value 
for which separability held, then all higher values of /i would also yield separability. But 
this turned out not to be invariably the case. So, it appeared that we needed to check the 
criterion 2,000 times for each point.) 



IV. CONCLUDING REMARKS 

In our earlier study we had also employed the Bloore parameterization of the two- 
qubit (and qubit-qutrit) systems to study the Hilbert-Schmidt separability probabilities of 
specialized systems of less than full dimensionality. We also reported an effort to determine 
a certain t/iree-dimensional function (in contrast to the one-dimensional functions f re ai{p) 
and fcompiex(n) above, but for somewhat a similar purpose) over the simplex of eigenvalues 
that would facilitate the calculation of the 15-dimensional volume of the two-qubit systems 
in terms of (monotone) metrics — such as the Bures, Kubo-Mori, Wigner-Yanase,. . . — 
other than the (non-monotone \4M ) Hilbert-Schmidt one considered here. (The Bloore 

, used above, did not seem immediately useful in this monotone metric 
context, since the eigenvalues of p are not explicitly expressed (cf. 44 1). Therefore, we had 
recourse in [^(J to the Euler-angle parameterization of Tilma, Byrd and Sudarshan 

Let us direct the reader to some papers of R. Kellerhals concerned, among other items, 
with the volumes of hyperbolic polyhedra (cf. 01 )■ In this line of work, the 

dilogarithm and, more generally, the polylogarithm functions play important roles. There 
have been some indicators in our investigations above (in particular, in integrations of the 
jacobians) that these functions also may be of relevance in our context. 

The extension to qubit-qutrit pairs (and even higher-dimensional compositie systems) of 
the univariate-function-strategy we have pursued above, for the case of qubit-qubit pairs, 
seems problematical, although we have not yet examined the matter in great detail. In 
the qubit-qubit case, the analysis is facilitated by the fact that it is sufficient that the 
determinant of the partial transpose be nonnegative for the Peres-Horodecki separability 
criterion to hold |36j, Thm. 5]. More requirements than this single one are needed in 
the qubit-qutrit scenario — even though the criterion on the nonnegativity of the partial 
transpose is still both necessary and sufficient for 6x6 density matrices. (In addition to the 
determinant, the leading minors and/or the individual eigenvalues of the partial transpose 
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of the 6x6 density matrix would need to be tested for nonnegativity, as well. Also the 
qubit-qutrit analogue of the ratio (fi) of diagonal entries would have to be defined, if even 
possible.) 
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